ON CONTINUOUS BRANCHES 
OF VERY SINGULAR SIMILARITY SOLUTIONS 
OF THE STABLE THIN FILM EQUATION 



J.D. EVANS AND V.A. GALAKTIONOV 

Abstract. The fourth-order thin film equation (the TFE-4) 

ut = -V • (M n VAu) + A(|u| p-1 u), where n > 0, p > 1, 

with the stable second-order diffusion term is considered. For the first critical exponent 

p = p Q = n + l + j ! for n £ (0, §), 

where N > 1 is the space dimension, the standard free-boundary problem (FBP) with 
zero height, zero contact angle, and zero-flux conditions is shown to admit continuous 
sets (branches) of source-type very singular self-similar similarity solutions (VSSs), 



u{x,t)=t 4 + ,iJV /(y), y = x/t< L +™ N . 

For the Cauchy problem (CP), continuous branches of oscillatory self-similar patterns of 
changing sign, which become "limits" of countable sets of FBP ones, are identified. 

For p ^ po, the set of VSSs is shown to be finite and to consist of a countable family 
of p-branches of similarity profiles that originate at a sequence of critical exponents 
{pu I > 0}. At p = pi, these p-branches appear via a nonlinear bifurcation mechanism 
from a countable set of similarity solutions of the second kind of the pure TFE 

ut = -V • (|u| n VAu) in l"xR + . 

Such solutions are detected by a combination of linear and nonlinear "Hermitian spec- 
tral theory" (in both the CP and FBP settings), which allows us to apply an analytical 
n-branching approach. This means constructing a continuous path as n — >• + to eigen- 
functions of a linear rescaled operator for n = 0, i.e., for the bi-harmonic equation 
Ut = — A 2 u. Numerics are used to confirm several analytical conclusions, which do not 
admit a fully rigorous study. 



1. Introduction: the stable TFE and main results 

1.1. The model and preliminary discussion. This paper is devoted to the study of 
large time behaviour of solutions of higher-order quasilinear degenerate parabolic equa- 
tions of not-divergent forms. More precisely, we construct global in time self-similar 
very singular solutions (VSSs) of the fourth-order quasilinear parabolic thin film equation 
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(TFE-4) with the stable homogeneous second-order diffusion term, 

(1.1) u t = -V ■ 0| n VAu) + A(|w| p-1 «), where n > and p > 1. 

The present results complete the analysis of the TFEs performed in [HI HO] , where count- 
able sets and continuous branches of blow-up and global similarity solutions were obtained 
for the limit unstable TFE with the backward diffusion parabolic term, 

(1.2) u t = -V • (\u\ n VAu) - A(\u\ p ^u) (n>0, p > 1). 

The main mathematical approaches for (11.21) also exhibit certain similarities with those 
applied in [TTJ to the sixth-order limit unstable TFE-6 

(1.3) u t = V ■(\u\ n VA 2 u)-A(\u\ p - 1 u) (n>0, p>l). 

Surveys and extended lists of related references on the physics and mathematics of such 
thin film PDEs can be found in these papers [TT] and [9] • As usual, for our future analysis, 
both pioneering papers of higher-order nonlinear diffusion theory in the 1990s by Bernis- 
Friedman [2] (mainly, FBP theory for TFEs), Bernis pQ and Bernis-McLeod jl], where 
oscillatory similarity solutions of the Cauchy problem for the fourth-order porous medium- 
like equations (the PME-4) were constructed, are key; see the monograph [261 Ch. 4] for 
further details in these areas. 

We consider for (11. ip the standard FBP with zero-height, zero contact angle, and con- 
servation of mass (zero flux) conditions, 

(1.4) u = Vu = -n ■ (\u\ n VAu - V(|u| p_1 u)) = 

at the singularity surface (interface) To [u] , which is the lateral boundary of supp u C 
M. N x M + with the unit outward normal n. For the range of n 6 (0, |) plus the semilinear 
case n = also discussed here, these three conditions are expected to give a correctly spec- 
ified problem for the fourth-order parabolic equation, which is completed with bounded, 
smooth, integrable, compactly supported initial data 

(1.5) u(x, 0) = u (x) in T [u] n {t = 0}. 

We will also treat the Cauchy problem (CP) for (11. ip with compactly supported initial 
data fll.5p in M. N . The CP admitting oscillatory solutions of "maximal" regularity will 
need a special setting. 

We begin our study in the critical "conservative" case 

(1-6) p = p = n + l + ±, 

which is easier technically and reveals specific properties of similarity patterns. Eventually, 
we extend our approach to p ^ p (more precisely, for p < p ). Notice that, for n = 0, 
equation (11.11) is the limit stable Cahn-Hilliard equation 

(1.7) u t = -A 2 u + A(\u\ p ~ l u), 

which occurs in various applications; see references in [T2] . 
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1.2. Main results and layout of the paper. We construct very singular self-similar 
(or source- type for p = p ) solutions of ( II ,ip in certain ranges of the parameters n, p 
and N. For small enough n > 0, we will often refer to analogies with the semilinear 
Cahn-Hilliard equation (II. 7ft . Typically, we assume that 

(1.8) n G (0, f) and p>n + l, 

but sometimes we also treat n > |, for which the CP continues to admit sign changing 
solutions that are infinitely oscillatory at the interfaces. 

In Section [2] we formulate the similarity setting of the problem. Our conclusions and 
further layout of the paper are as follows: we show that the stable TFE (11.11) admits: 

(i) In the critical case p = po, continuous families of global similarity solutions of the 
CP (Section ED and of the FBP (Section EJ); 

(ii) As a co-product, we first study countable set of similarity solutions of the pure 
TFE (these define special bifurcation values {p{\ for the full model (II. ip . Section S]) 

(1.9) u t = -V ■ (|u| n VAu) in M N x R + - and 

(iii) For p < po (for p > p no such VSSs exist), the number of similarity solutions (for a 
given p value) becomes finite, and in the CP, there exists a countable family of p-branches 
of similarity profiles that originate at certain nonlinear bifurcation points {pi > 1, I > 0} 
(Section [5]). Some of the results are extended to the FBP (Section [9]). 

The principal issue that occurs is the actual relation between similarity solutions of 
the standard FBP and the Cauchy problem (the latter are infinitely oscillatory near the 
interfaces for not that large n). More clearly and convincingly than in our previous 
research, we show that, in several cases, for each asymptotic pattern of the CP, there 
exists a countable set of FBP patterns, which eventually converge to the CP one. It turns 
out that this is a rather general principal even for the linear problem for n = 0, i.e., for 
the bi-harmonic equation, 

(1.10) u t = —A 2 u. 

To show this, we will need to develop a type of Hermitian spectral theory for linear rescaled 
operators for both the CP (Section l4.ip and for the FBP setting (Section l9.2p . The latter 
one becomes more difficult and unusual, with a multi-dimensional space of eigenvalues, 
so we discuss just initial aspects of such a theory therein. 

Similar principles of self-similar asymptotics apply [TTJ to the sixth-order stable TFE 

(1.11) ut = V • (\u\ n VA 2 u) + A{\u\ p - l u), 

for n G [0, |). In this case, the first critical exponent is po = n + 1 + The semilinear 
case n = leads to the sixth-order unstable limit Cahn-Hilliard equation 

(1.12) Ut = A 3 m + A(|"u| p_1, u), 



whose similarity solutions can be studied as in [T2] . 
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2. Global similarity solutions: general statement and preliminaries 
The similarity solutions of (ll.ip have the form 

(2.1) u s (x, t) = r a f(y), y = x/t p , with a = > 0, = ^+2)] > °- 
The function / solves a quasilinear elliptic equation, namely, 

(2.2) A+(/) = -V • [|/| B VA/ - Vd/r 1 /)] + ft, • V/ + af = 0. 

For n > 0, a natural functional setting for both the FBP and the CP includes the condition 

(2.3) f(y) is non-trivial in a bounded domain in M. N . 

In the CP, f(y) can be extended by f(y) = outside the support. For FBPs, posed 
by definition in a bounded domain, such an extension is not applicable. Note that, for 
the CP, the elliptic equation ( 12.21) admits non-compactly supported solutions with the 
asymptotics as y — > oo governed by the leading linear first-order operator: 

(2.4) /3y-V/ + a/ + ... = f(y) = C\y\-f(l + o(l)), 

where C = C(i^r) is an arbitrary smooth function on the unit sphere S 1 ^ -1 . Actually, in 
order to satisfy the desired condition (12. 3p . one needs to demand that 

(2.5) C = in (JMD- 
For the CP with n = 0, is replaced by (see [12] ) 

(2.6) /(y) has exponential decay at infinity, 

meaning that / belongs to a special weighted L 2 space. The condition (12. 5p is also 
necessarily implied. 

Under (12. 3p (or equivalently (12. 6p ). integrating (12. ip over M. N yields the following mass 
time-dependence for p ^ p : 

n (p-pq) 

(2.7) fus(x,t)dx = t*&=&toff(y)dy jf = {p^ Po ). 

For p = po, any mass of f(y) is formally allowed; cf. [ 10] . 

For general solutions of the TFE-4 fll.ip . the self-similar scaling 

(2.8) u(x,t) = (l + t)~ a 6(y,r), y = x/(l + t)^ r = ln(l + t), 
yields the evolution equation with the same operator as in (12.21) . 

(2.9) 6 T = A + (6) for r > 0. 

Then a typical asymptotic stabilization problem occurs asr-y +oo, which in particular, 
requires to know all possible equilibria of A + , this being the main current problem of 
concern. 

The critical exponent (II .6p follows from conservation of mass; q.v. the same derivation 
in [91 § 3]. In addition, a countable sequence of other critical exponents {pi, I = 0, 1, 2, ...} 
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is expected to exist. This is confirmed in the semilinear case n = (for the limit Cahn- 
Hilliard equation (11.7)) ) . where [121 § 5] 

(2.10) Pi = i + _2_ f or any / = 0,1,2,.... 

See further comments in [9| § 2]. 

3. The CP: local oscillatory ODE bundles and profiles for p = p 

We next study the similarity ODE in the radial setting. Let y denote the single spatial 
variable \y\ > 0. The operator of the equation (12.11) is then ordinary differential, 

(3.1) A+(/) = [v N - 1 \f\ n (^{v N - 1 n , y ~ y N - 1 (\fr 1 f)'} , + 0vf + af = 0. 

We first describe the corresponding oscillatory bundle of asymptotic profiles close to 
interface points, which are attributed to the CP. It turns out that we have to begin with 
the CP. Namely, we will show that similar FBP profiles are naturally associated with the 
CP ones and, in particular, for both n = and n > demand certain extra "Hermitian 
spectral theory" as an extended and more difficult version of that in the Cauchy setting 
in R N . 

3.1. The Cauchy problem: local oscillatory behaviour close to interfaces. These 
questions have been considered before, so, following [10, § 7.1], we briefly indicate the 
oscillatory asymptotic bundle of similarity profiles f(y) exhibiting maximal regularity at 
the interface y = y , so that, being extended by / = for y > y , these will give solutions 
of the CP. It is easy to see that, for n e (0, |), the ODE (13.21) does not admit nonnegative 
solutions of the maximal regularity, which can be considered as a counterpart of smooth 
similarity solutions of the CP for n = [121 § 5], i.e., admitting a regular limit as n — > + . 
Hence, the ODE ( 13.2)) implies that sufficiently regular solutions f(y) must be oscillatory 
near interfaces; cf. proofs in [3]. It is important to prescribe the precise structure of 
such oscillatory singularities of the ODE and determine the dimension of the asymptotic 
bundle. 

For n > 0, we take the thin film ODE (13. ip . keeping the main terms for y fa y^ and 
integrating once, to obtain 

(3.2) |/|« (/" + )' - (l/r 1 /)' = 0Vof + (higher-order terms) . 

For N = 1, choosing next just two leading terms close to the interface yields (in fact, one 
can see that this approximation holds for any dimension N > 1) 

(3.3) | / |« / ,,/ = /3i// + ... = Ao/ + ..., 
where Ao = (3yo- Thus, we need to consider the unperturbed ODE 

(3.4) 1/ ! n f = Ao/ (A = /3y >0), 

whose orbits will be exponentially small perturbations near interfaces of that for ( 13.11) . 
The ODE (13. 4p has the following representation of the solutions [TUl § 7]: as y — > t/q", 

(3.5) f(y) = (y - yY<p(v) , V = Mvo~y), with fi = |, 
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where the oscillatory component <p satisfies the autonomous ODE 

(3.6) iff" + 3(/i - l)<p" + (3u 2 - 6u + 2)ip' + n{n - 1)0* - 2)<p + A M~> = 0. 

One can see that, on orbits like (I3.5p . the term /)', which has been neglected in 

(13. 2p . is much small than / for all p > 1 + | that is true since always p > 1 + n. 

We are interested in periodic solutions of (I3.6p . which will give, according to (|3.5|) . 
oscillatory profiles changing sign infinitely many times as y — > y^, i.e., as rj — > — oo. 
Via (|3.5|) . periodic functions v ? *( ? ?) establish the simplest oscillatory connections with the 
interface points keeping the maximal regularity of the envelope: 

3 

f(y) ~ (vo-y) 5 for y^y^, 

which represents the true scaling- invariant nature of the ODE (13. 4p . According to (13. 5p . 
the regularity at y = yo then increases as n — > + forming, at n = 0, analytic solutions. In 
[TTj § 6] , we present a discussion and references related to the theory of periodic solutions 
of higher-order ODEs. Unlike the fifth-order case in [TTI § 4] , the ODE (13. 6p is of third 
order and can be reduced to a first-order ODE; see [TBI § 7.1]. Therefore the existence 
of a periodic solution is not principally difficult, while the uniqueness (and the stability) 
are harder. We expect, and this is confirmed by numerics [TO], that this limit cycle is 
"almost" globally stable (note that all the orbits of (13.61) are uniformly bounded, so a 
stable attractor should be available, though sometimes may have a stable manifold, 
which was not observed) and is unique. 

Increasing n more, this periodic solution ^p*(s) is destroyed in a heteroclinic bifurcation 
at the point (XQJ § 7.2] 

(3.7) n h = 1.758665... (and n h e (|,n+), where n + = ^= = 1.9019238... , ^7\), 

with a standard scenario of homoclinic/heteroclinic bifurcations, [2U Ch. 4]. A rigorous 
justification of such non-local bifurcations is an open problem. 

Thus, for n larger than |, not all the solutions are oscillatory near the interfaces. 
For n G (|,3), there exists a one-parametric bundle of positive solutions with constant 
equilibria <p(rj) = ±ip given by 

(3-8) ±<p = ±[- ^Jt>l_ 2) ]K 

For matching purposes, this is not enough and the whole 2D asymptotic bundle (13. 5p of 
oscillatory solutions has to be taken into account, so, for the CP, the oscillatory behaviour 
is expected to remain generic for all n e (0, nh) (and similar to the linear case n = with 
the interface at y — oo). 

Thus [TU], it is key that, for n G (0, h^), there exists a 2D bundle of asymptotic oscilla- 
tory orbits near the interface with the behaviour (here yo > and sq G R are parameters) 

(3.9) f(y) = (y -y)n( P:> ,(ln(y -y) + s ) + ... as y -> y . 

Another important question is the passage to the limit n — > + that shows convergence 
to solutions of the semilinear Cahn-Hilliard equation. This is explained in detail in [T0~| 
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The Cauchy problem: N=1 , p=p =n+3, n=1 ; parameterization by f(0) 



3 
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f(y) 
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1 

0.5 
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y 

FIGURE 1. Similarity profiles for the CP as solutions of (f37l0|) . (pTTTj) for N = 
n = 1, p = n + 3 = 4. 

§ 7.6]. Various oscillatory sign change issues for nonlinear degenerate higher-order PDEs 
of different types are addressed in [181 Ch. 3-5]. 

3.2. Continuous branches of similarity profiles for p = pq. Again, without loss of 
generality, we treat the case N — 1, where the ODE is simpler and takes the form 

(s.io) \f\ n f"'-^yf-(\f\ n fy = o. 

The origin of the existence of continuous branches (parameterized by e.g. mass) is the 
fact that the ODE (I3.10p is of third order. Thus the single symmetry condition 

(3.11) f(0)=0 (/(0)^0). 

is posed at the origin, while the shooting bundle from the singularity point y = yo is 2D 
according to (13. 9p . So one parameter is free. The existence by a shooting approach is 
standard, so we refer to [TUJ [12[ [19] as a guide to such equations. 

The numerical results below were mainly obtained using Matlab's two-point boundary 
value problem collocation solver bvp4c (with default parameter values RelTol= 10 -3 , 
AbsTol= 1CT 6 ). The standard regularization 

(3.12) |/|« (5 2 + f 2 )" 

was used with 5 = 1CT 6 (and sometimes up to 10~ 10 in order to see complicated and 
refined zero structure of solutions) for q = n and j?. 

In Figure [H we present similarity profile for n = N = 1 and p = po = n + 3 = A, which 
are parameterized by the values at the origin /(0). 
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The Cauchy problem: N=1 , p=p =3, n=0; parameterization by f(0) The Cauchy prob|em: comparison of 2er0 sets? N=1 n=1 and n=0 




5 10 15 20 25 ' 10 15 20 25 

y y 

(a) f(y) for n = (b) zero structure enlarged 

FIGURE 2. Similarity profiles for the CP satisfying (j330|) . (f37TT]) for iV = 1, 
n = p = 3; profiles (a) and zero structure (b). 

For comparison, in Figure E^a), we present similar profiles for the semilinear case n = 0, 
i.e., for the limit CH equation fll.Tj) . where p = po = 3. FigureE^b) shows a clear difference 
of the "tail" behaviour for n = 1 (nonlinear oscillations (13. 9j) ) and n = (a linearized 
behaviour) . 

In Figure [3j we explain how the similarity profiles at p = po = n + 3 are deformed with 
n starting again from the semilinear case n = 0. It is seen that the profiles get thinner as 
n increases and also get less oscillatory near the interface. Note that the last case with 
the maximal n = 1.75 (the CP profile must still change sign, while the FBP one does not 
[TO] ) is close to the heteroclinic value (13. 7p . after which the similarity profiles are assumed 
to be finitely oscillatory, i.e., can have a finite number of sign changes near the interface 
(or no at all); see further comments in [10, § 9.4]. This finite oscillation phenomena lead 
to difficult and challenging numerical problems. Our numerics show that, for n = 1.75, 
the CP profile still changes sign near the interface; see the next Figure HI 

In Figure |3l we also included the case of a single negative n = — |, which gives a standard 
source-type profile but, of course, without finite interface, where f(y) is oscillatory as 
y — > +oo. The structure of those oscillations at infinity is difficult and different from 
already known ones, and is not studied here. 

4. The Cauchy problem: countable family of source-type profiles via an 

n-branching approach 

For future convenience, we now need to postpone our study of the original PDE (11.11) 
and digress to the pure unperturbed TFE. In general, construction of various oscillatory 
source-type solutions of the CP for pure TFE (11.91) is a difficult nonlinear problem, which 
is harder than that for the FBP. 
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The Cauchy problem: N=1 , p=n+3, n=1 .75, 1 .5, 1 , 0.5, 0, -0.5 




FIGURE 3. Similarity profiles for the CP as solutions of (|57T0|) . (ETTT]) for N = 1, 
p = n + 3 = 4 and various n G [—5, if]; parameterisation is /(0) = 1. 



The Cauchy problem: N=1 , p=n+3, n=1 .75, 1 .5, 1 , 0.5, 0, -0.5; enlarged 




FIGURE 4. Enlarged zero structure of VSS profiles from Figure [3] for p = 1.75, 
f(y) still changes sign. 

4.1. The linear case n = 0: basics of Hermitian spectral theory. On the other 
hand, for n = 0, i.e., for the bi-harmonic equation f ll.lOp . the first such profile exists, 
is unique (up to mass scaling), and is just the rescaled kernel F(y) of the fundamental 
solution of (JTIDJ) : 



(4.i; 



BF 



b(x,t)—t *F(y), y = x/t-i, where 
-A 2 F+±y- VF + f F = in R N , J F = 1; 



see 0, § 4]. Moreover, there exists a countable set of eigenf unctions {ip 7 , I — \j\ > 0} of 
the corresponding rescaled non-selfO-adjoint operator B with the discrete spectrum [8] 

(4.2) a ( B ) = {-I, 1 = 0,1,2,...}. 
The eigenfunctions are derivatives of the rescaled kernel F, 

(4.3) f(t,)^^D^) for any multiindex 7. 
In addition, the adjoint operator has the same spectrum: 

(4.4) B * = -A 2 -\yV with a(B*) = a(B) = { - I I = 0, 1, 2, ...}, 

and a complete set of eigenfunctions {^*(?/)}, which are generalized Hermite polynomials. 
See more details on such a Hermitian spectral theory of the pair {B, B*} in [HI | 



4.2. n-branching of similarity solutions. We apply the n-branching approach from 
the linear case n = in order to explain existence of a countable set of similarity solutions 
of the TFE (11.91) . We will follow classic branching theory in the case of non-analytic 
nonlinearities of finite regularity; see [251 § 27] and [231 Ch. 8]. 

We look for solutions of (11.91) with small n > in the standard form 

(4.5) u,(x,t)=t- a f(y), y = x/tP, where = 

where the multiindex 7 is used for numbering to be explained later on (in fact, similar to 
the linear eigenfunctions (14. 3j) ) . Then / = / 7 (y) solves the elliptic equation 

(4.6) A„(/) = -V.(|/| n VA/)+/fy.V/ + a/ = in R N . 
Note that, in general, for I > 1, we have to assume that 

(4.7) Jf(y)dy = 0, 

so that the solutions ( 14.5ft satisfy the mass conservation condition 

(4.8) fu 1 (x,t)dx = (|7| > 1). 

For / = 0, where a = (3N and (3 = 4 + nN , the assumption (14.71) is not necessary, since the 
PDE (14. 6 p is fully divergent and admits integration over M. N . 
For small n > in ( 14.6p . we have 

(4.9) /9 = i-fn, 
and use the following expansion: 

(4.10) |/r = 1+nln 1/1 + o(n). 

Here, ( I4.10p should be understood in the weak sense, which is necessary for using in the 
equivalent integral equation; see Proposition 14.11 below. 

Substituting expansions (14. 9p and (I4.10p (still completely formal) into ( 14. 6 p yields 

(4.11) A n (/) = B/ + (a - f )/ + nC(f) + o(n) = 0, 
with the perturbation operator 

(4.12) £(/) = -V • (In l/l VA/) - f j/ • V/. 
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We next describe the behaviour of solutions for small n > and apply the classical 
Lyapunov-Schmidt method [231 Ch. 8] to equation (14.111) . Recall that, in this linearized 
setting, we naturally arrive at the functional framework that is suitable for the linear 
operator B, i.e., it is lA(R ), with the domain H*(M. ), etc., and a similar setting for the 
adjoint operator B*; see above and further details in [8]. 

Therefore, for n = 0, we have to study branching of a nonlinear eigenfunction from the 
linear one, where / is a certain nontrivial finite linear combination of eigenfunctions from 
a given eigenspace with fixed A 7 = — = — |, i.e., 

(4-13) f = 0i = Z hH C^ (^0). 

Those eigenfunctions are just derivatives (14.3p of the analytic radially symmetric rescaled 
kernel P(|y|) of the fundamental solution ( 14. ip . Therefore, the nodal (zero) set of such / 
in (14.131) is well understood and consists of a countable set of isolated sufficiently smooth 
hypersurfaces which can concentrate as y — > oo, where 

(4.14) 4>i{y) as y — > oo uniformly and exponentially fast. 
Hence, returning to the key limit (I4.10p . we have: 

Proposition 4.1. For a function f given by l \4- 15fy , in the sense of distributions, 

(4.15) i (l/l" -1)- In |/| as 0+ 

According to ( 14.15p . analyzing the integral equation for /, we can use the fact that, for 
any function G C (and G Co(^ N )), 

(4.16) f R A\f(y)\ n - dy = n[J RN In |/(y) dy + o(l)]. 

It follows from (14. lip that branching is possible under the following non-trivial kernel 
assumption: for n = 0, 

(4.17) a -N = _ k= i_ ^ oi(0) = ^, l>0. 

This gives the countable sequence of critical exponents {ai(n), (3i(n)} (to be determined) 
of the similarity patterns (14.51) of the TFE for small n > 0. 
By [S] Lemma 4.1], the kernel of the linearized operator 

E = ker (B - A,/) = Span {^, \/3\ = 1} 

is finite-dimensional. Hence, denoting by E\ the complementary (orthogonal to ^o) in- 
variant subspace, we set 

(4.18) / = 0/ + Vi, where (pi G £ and V 1 = $} 7 [>j c 7^7 E E i- 

According to the known spectral properties of operator B, we define and Pi, Po + Pi — 
I, to be projections onto E Q and E\ respectively. We also introduce a perturbation of the 
parameter a by setting 



(4.19) 



oti(n) = ai(0) + 5, with 5 = 5(n). 
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This perturbation 5 is obtained from the orthogonality condition by substituting into 
(14. lip and multiplying by ip* This gives 

(4.20) 5(n) = c t n + o(n), 
where q is obtained from the system 

(4.21) (/m^;> = c h h\ = L 

Since according to (I4.18P , (pi is given by (14.131) , (I4.2ip is an algebraic system for unknowns 
{C 7 } and q. It can be solved, for instance, in the radial geometry and in some other cases 
(including those where the dimension of the kernel is odd; even dimensions are known to 
need additional treatment); see more details in [T6| App. A]. However, the total number 
of solutions of the non- variational system (14.211) remains unclear. 
Finally, setting 

(4.22) V x =nY + o(n), 

we obtain, passing to the limit n — > + , the following equation for Y: 

(4.23) BY = -c l <p l +C(<j> l ). 

By Fredholm's theory, in view of the orthogonality, it admits a unique solution Y £ E\. 

In general, the above analysis shows that, up to solvability of the nonlinear algebraic 
systems, the TFE admits a countable set of different source-type similarity solutions ( 14 . 5 p 
at least for small n > 0, where the parameters ai(n) are given by 

(4.24) ai{n) = ^ + cin + o{n) as n ->• 0+; 1 = 0,1,2,.... 

At n = 0, these solutions are originated from suitable eigenfunctions of the linear operator 
in ( 14.11) . The global extensions of these n-branches of similarity solutions for larger n > 
represent a difficult open problem, to be treated numerically later on. 

4.3. Nonlinear eigenfunctions of the TFE in one dimension. We consider the 
Cauchy problem for the ID TFE with continuous compactly supported initial data, 

(4.25) u t = -{\u\ n u xxx ) x in R x R + , u(x, 0) = u {x) G C (R). 

Then, for A = 1, the nonlinear eigenvalue problem for the elliptic equation (14.61) is 
formulated as follows: 



(4.26) -(\f\n f »y + l=«ny f > + af = in M> f( y )^o, f £ Co 



Actually, (I4.26P is also about self-similarity of second kind, where the desired set of 
parameters (nonlinear eigenvalues) {ai(n), I > 0} is obtained not by a pure dimensional 
analysis, but via solvability of a nonlinear ODE in a given functional class C*o(R) of 
compactly supported functions satisfying the condition of maximal regularity. The term 
similarity of the second type was introduced by Ya.B. Zel'dovich in 1956, [2T] . 

Note that, for n = 0, (I4.26P in L 2 (R), where we replace the last condition by 




is a standard linear eigenvalue problem for a non self- adjoint operator with the point 
spectrum ( 14. 2 p and complete-closed set of eigenfunctions given in ( 14.31) . jH]- 

The first nonlinear eigenvalue-eigenf unction pair {Fo, «o} of ( 14.26}) was proved to exist 
for n G (0, 1]; see [TOj § 9]. In this case, the first eigenvalue is 

(4-27) «oH = 4T^L=i = 4T^ 

It turns out that, in view of the highly oscillatory nature of those profiles near interfaces, 
even identifying the position of interfaces numerically, is not an easy problem. Therefore, 
we begin with Figure |5j where the first even nonlinear eigenfunction is presented for n = 0, 
~, and 1. A careful study of their zero structure in the log-scale in (b) allows us to find an 
approximate and rather rough interface location, according to the expansion (13. 5p . which 
yields 

(4-28) In \f(y)\ = f ln(y -y) + ln \<p(Hy -y)\ + .... 

For n — 0, the expansion is exponential (cf. ( 19 .9p below) and is entirely different, which 
is seen in (b); recall the regularization such as in (I3.12p eventually entering the expansion 
for \ f\ very small. 

Thus, in what follows, we use the parameterisation: 

(4.29) F,(0) = 1, 1 = 0,2,4,...; F/(0) = 1, Z = l,3,.... 
We then obtain that, respectively, 

(4.30) a x = 0.2534..., a 2 = 0.320... (n=l). 

In addition, numerics show that «3 ~ 0.38 for n = 1. 

The results of an accurate numerical study of the first four nonlinear eigenfunctions for 
N = 1 and various n e [0, 2] are presented in Figure [61 Further eigenfunctions are very 
difficult to obtain numerically, to say nothing about an analytical proof of their existence. 
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In Figured we show a formal schematic behaviour of the nonlinear eigenvalues cfy(n) 
of (I4.26p . The first n-branch, according to (I4.27p . has the explicit form 

®o(n) = jjjj, n G [0,3). 

Other n-branches in Figure [7J are not explicit and are hypothetical. According to the re- 
branching approach, all these branches originate at the eigenvalues of the linear problem 
(O, i.e., 

(4.31) a,(0) = -A J+1 = for Z = 0,l,2,..., 

and moreover, after scaling, we may assume that, at n — 0, the similarity profiles Fi(y) 
coincide with the eigenfunctions ( 14.3}) . and hence by continuity mimic their geometric 
shapes for n > 0. 

Figure El shows the actual numerical construction of first four re-branches, and even 
these involve technical difficulties. Note that, at the critical heteroclinic bifurcation value 
(13. 7p . the similarity profiles Fi(y) are supposed to loose their oscillatory behaviour at the 
interface and become finite oscillatory for n > rih (or even non-oscillatory at all); see [T0| 
§7.2]. 

Analytical difficulties for the eigenvalue problem (I4.26P begin already with I = 1, i.e., 
with the dipole profile Fi(y). This study has a well-developed history (see [3j El [7J and 
references therein), but still there are no definite results of existence and uniqueness of 
F\ in both the FBP and Cauchy problem settings. 

We end this discussion with the following: 

Conjecture [4J.1. (i) For any n e (0, n^), the nonlinear eigenvalue problem { 4-26 ) admits 
a countable set of sufficiently smooth solutions of maximal regularity 

(4.32) ® = {F l (y), 1 = 0,1,2,...}, 

where the nonlinear eigenvalues {a{\ form a strictly increasing sequence and 

(4.33) ai — > - as I — > oo. 

(ii) The eigenfunction subset { 4-32} is evolutionary complete in Cq{M) for the TFE {4-25), 
i.e., for any u ^ 0, there exists a finite I > and a constant b = b(u ) ^ such that 

(4.34) u(x, t) = r ai [bFix/t? 1 \b\ n/i ) + o(l)] as t -)• oo. 



We expect that an analogous countable set of radially symmetric similarity solutions 
exists for the TFE (11.91) in any dimension N > 2, though numerical calculations become 
much more difficult than for N = 1. Moreover, the branching approach in Section I4T21 
shows that there many other non-radial similarity solutions that have a more complicated 
geometry but which for small n > mimic the eigenfunctions (14. 3p . 

The evolution completeness of nonlinear eigenfunctions is known rigorously for the PME 
(4.35) u t = {\u\ n u) xx (n>0) 



More details on this are given in [TO] . 
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The second (dipole) profile F (y), N=1 



The first (ZKB) profile F (y), N=1 




The third profile F (y), N=1 





The fourth profile F (y), N=1 











n=0 






- - n=0.5 






n=1 






- - n=1.5 






n=2 










V /■ '• \ / // 



FIGURE 6. Illustrative numerical solutions for the first four nonlinear eigen- 
functions -F}(y), I = 0, 1, 2, 3, shown for selected n in one dimension N = 1. The 
corresponding behaviour of the eigenvalues are shown in Figure [HJ The regular- 
isation (|3.12p with 5 = 10~ 2 was used in the numerical shooting scheme, where 
the even numbered profiles satisfy i*}(0) = l,F/(0) = F["(0) = 0, whilst the odd 
numbered profiles have F/(0) = 1, F,(0) = F/'(0) = 0. 



in a bounded interval [14]; see also [16] for results in M. N for initial data uq G Cq(R. n ) and 
on an n-branching technique. Existence of a countable set of radial similarity solutions of 
(Q5]l in R N x R + was proved by Hulshof 



5. The Cauchy problem: on nonlinear ^-bifurcations 

5.1. Semilinear Cahn— Hilliard equation: countable set of critical exponents. In 

order to explain the essence of the nonlinear bifurcation analysis, we first digress to the 
stable CH equation (11. 7p . for which the analysis is much simpler; cf. [121 EH]- 
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nonlinear eigenvalues a (n) for F(y), 1=0,1,2,... 




FIGURE 7. Formal n-branches of nonlinear eigenvalues ai(n) of (I4.26D . 



Nonlinear eigenvalues cc^n) for F^y), 1=0,1 ,2,3 




FIGURE 8. The actual n-branches of the first four nonlinear eigenvalues ai(n), 
I = 0,1, 2, 3, constructed numerically. 

Namely, studying the behaviour as t — » +oo, we perform the standard scaling 
(5.1) u(x,t) = (1 +ty^=^v(y,r), y = x/(l + t)* } r = m(l + t), 



where v(y, r) solves the following rescaled equation: 

(5.2) v T = -A 2 v + \y-Vv + v + A^v^v) = Bv + c v + AQv^v). 

It then follows from (14. 2 p that a centre manifold behaviour is formally possible in the 
critical cases (12.101) only: 

(5-3) Co=2(FiJ-f = l P = Pi = l + jki (^°)' 

Checking the necessary condition of such a centre manifold behaviour and looking, say, 
for a solution moving along the centre eigenspace, 

(5.4) v(y, t) = a 7 (r)^ 7 (j/) + w, W-Lip 7 , sup y \w(y, r) | = o(a 7 (r)) as r — )• oo, 
and substituting into (15. 2p yields on multiplication by ip* (see [8] for details) 

(5.5) d 7 (r) = /i 7 |a 7 | p_1 a 7 + where /i 7 = (A\ip 1 \ p ~ 1 ip 1 , ip*) = (|'0 7 | p ~ 

Since ^-(y) is a 7-degree polynomial [8], we then conclude that the necessary condition 
of existence of such a centre subspace behaviour is as follows: 

(5.6) yU 7 7^ at least, for I7I > 2. 
Note that, in the limit p — > 1, the following holds: 

(5.7) /x 7 = 1 by bi-orthonormality of eigenfunctions, 

so that (15. 6 p is true for / ^> 1 by continuity of the integral relative to the parameter p. 

Eventually, the centre subspace behaviour (I5.5P generates the following asymptotic 
patterns for the CH (fTTjl : 

(5.8) u 7 (x,t) ~ C 7 (tln 2 fr^^ 7 (^) + ... as t -> 00 (Z = \-y\ > 2), 
where constants C 7 are independent of initial data Uq. 

5.2. Local bifurcations from p\. We now return to the VSSs of the TFE with the stable 
PME term (II. ip and perform a formal nonlinear version of a ^-bifurcation (branching) 
analysis for n > 0. As usual, according to classic branching theory [231 122], a justifica- 
tion (if any) is performed for the equivalent quasilinear integral equation with compact 
operators. For simplicity, we present computations for the differential setting. 

Thus, we consider the elliptic PDE (12. 2p . The critical exponents {p{\ are then deter- 
mined from the equality (q.v. ( 15.31) ) 

M) a ee 2pi[n) i in+2) = a*(n) ft(n) = *±2 + ^ (I > 0). 

In particular, for the semilinear case n = 0, we have a/(0) = ^ from (14TT7) . so that 
(15. 9p leads to (15.31) . i.e., to the known sequence of critical exponents (12.101) . 
We next use an expansion relative to the small parameter e = po — p, i.e., as e — > 0, 

a = 2w _ (n 1 +2 )-2 £ = a i + 2a2 i e + - > 
/3 = ^P + Q. + ..., ^^[n + 2 + i-^]. 
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Substituting these expansions and the last one in (14.91) into (12. 2p and performing the same 
standard linearization yields 

(5.10) A n (f) + A(|/|»/) + e[ - A(|/|»/ln |/|) + &F] + 0(e 2 ) = 0, 

where C\ — C\y - V + 2a 2 / 

is a linear operator, and A n is the rescaled operator (14. 6 j) of the pure TFE with the 
parameter a = a/(n) (an eigenvalue), for which there exists the corresponding similarity 
profile F\(y) (the nonlinear eigenfunction). The fact that the operator A n with a = ai 
in (15. lOf) occurs in the rescaled pure TFE correctly describes the essence of a "nonlinear 
bifurcation phenomenon" to be revealed. 
To this end, we use the additional invariant scaling of the operator A„ by setting 

(5.11) f(y) = bF(y/b%) (b > 0), 
where b = b(e) > is a small parameter satisfying 

(5.12) b(e) ^ as e 0, 

to be determined. Substituting (15.111) into (I5.10p and omitting all higher-order terms 
(including the one with the logarithmic multiplier In |&(e)|) yields 

(5.13) A„(F) + 6P J -f A(\F\ Pl F) + eC x F = 0. 

Finally, we perform linearization about the nonlinear eigenfunction Fi(y) by setting 

F — Fi + Y. 

This yields the following linear non-homogeneous problem: 

(5.14) K( F i) Y + 6»«-f A(|F,|«F,) + eA^i = 0. 
Here the derivative is given by 

A' n (F)Y = -V • (VAF)F + VAF)] + fry ■ VY + a^F. 

The rest of the analysis depends on assumed good spectral properties of the linearised 
operator A! n {Fj). We follow the lines of a similar analysis performed for the FBP case 
in [HJ § 2], where the operator A' n (Fi) for n = 1 turns out to possess a (Friedrichs') 
self-adjoint extension with compact resolvent and discrete spectrum. Such a self-adjoint 
extension does not exist for the oscillatory F(y). Here we use general theory of non- 
self-adjoint operators; see e.g., [20J. A proper functional setting of this operator is more 
straightforward for N — 1 (and in the radial setting), where, using the behaviour of 
F(y) — > as y — > 1, it is possible to check whether the resolvent is compact in a suitable 
weighted L 2 space. In general, this is a difficult problem; see below. 

We assume that such a proper functional setting is available for A n , so we deal with 
operators having solutions with "minimal" singularities at the boundary of the support 
Si, where the operator is degenerate and singular. Namely, we assume that A' n (Fi) has 
discrete spectrum and a complete and closed set of eigenfunctions denoted again by {^ 7 }. 
We also assume that the kernel is finite dimensional and we are able to determine the 
spectrum, eigenfunctions {i/jZ}, and the kernel of the adjoint operator (A' n (Fi))* defined 
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in a natural way using the topology of the dual space L 2 and having the same point 
spectrum (the latter is true for compact operators in a suitable space [221 Ch. 4]). 

Further, we assume that there exists the orthogonal subspace Span{^ 7 , 17I > /} of 
eigenfunctions of A' n (Fi), and we look for solutions of (I5.14j) in the form 

Y = <f>i + w, 

where <pi belongs to the kernel and hence is analogously given by (14. 13ft and w belongs to 
the orthogonal complement of the kernel. In doing so, we need to transform ( 15. 14j) into 
an equivalent integral equation with compact operators, but for convenience, we continue 
our computations using the differential version; see additional details in [19 , § 3]. 

Thus, multiplying (I5.14p by if)* with any I7I = I in L 2 and, if necessary, integrating 
by parts in the differential term y ■ Vi 7 ] in C1F1, we obtain the following orthogonality 
condition of solvability (Lyapunov-Schmidt's branching equation [251 § 27]): 

(5.15) ^-i(A(|F^- 1 F i ),V;> = -s(C 2 F h ^;) for all | 7 | = I. 

These are algebraic equations for the expansion coefficients {C 7 } in (I4.13p and the function 
b = b(e). Similar to (15. 5p . one needs to check whether the constants are non zero, 

(5.16) (Adi^- 1 ^),^;) ^0 and (AF^^O, 

which is not a simple problem and can lead to restrictions for such a behaviour. The 
analysis is much simpler if the kernel is ID, which always happens in the radial geometry 
where we deal with ordinary differential operators. Then (I5.15P is a single and easily 
solved algebraic equation, for which the "transversality" problem (I5.16P also occurs. 
Under the conditions (I5.16p . the parameter b(e) in (15. lip for p rj p is given by 

(5.17) b{e) ~ [ 7z (pj -p)}^. 

The direction of each p^-branch and, whether the bifurcation is sub- or supercritical, 
depends on the sign on the coefficient 7; that follows from (I5.15p . This can be checked 
numerically only, but, in general, we expect that the most of these nonlinear bifurcations 
are sub critical so the ^-branches exist for p < pi. 

For n = 0, a rigorous justification of this bifurcation analysis can be found in [191 § 6], 
where a countable number of p-branches was shown to originate at bifurcation points 
(12.101) and were detected on the basis of known spectral properties of the corresponding 
linear operator in (j4.ip ; see details in [8]. For n > 0, as we have seen, the justification 
needs spectral properties of the linearised operator A' n (Fi) and the corresponding adjoint 
one (A' n (Fi))*, which is very difficult for non-radial nonlinear eigenfunctions Fi and is 
an open problem. In particular, it would be important to know that the bi-orthonormal 
eigenfunction subset {^ 7 } of the operator A' n (Fi) is complete and closed in a weighted 
L 2 -space or in some specially defined closed subspace (for n = 0, such results are available 
[8]). We expect that for n«0, there exist critical exponents for the TFE with absorption 
that are close to those in (12.101) at n — 0. This can be checked by standard branching-type 
calculus; see [161 App. A], where nonlinear eigenfunctions of the rescaled PME in were 
studied by a branching approach. 
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6. The Cauchy problem: towards global extensions of ^-branches 

6.1. Examples of various profiles. Recall that, for p ^ p in the ODE ( 13. ip . we still 
have a 2D bundle at the singular interface point (13. 9p . but now, for even profiles, we also 
need to satisfy two symmetry boundary conditions at the origin: 

(6.1) f(0) = f" (0) = 0. 

Therefore, unlike the third-order problem ( I3.10p . ( 13. lip in the critical case p = po, we 
cannot expect continuous sets of solutions. Actually, as was shown in JTUJ [TTJ [121 UHL i n 
these non-critical cases, there occurs a countable set of p-branches of similarity profiles, 
which originate at the standard (for n = 0) or nonlinear bifurcation points {pi} as ex- 
plained in Section [5j The global behaviour of such p-branches can be complicated and we 
do not intend to study these delicate open questions in any detail, restricting ourselves 
to examples only. 

In Figure [HI we present some VSS profile f(y) for N = 1 in two cases: n = 1 and 
p = 3 < po = 4 in (a) and n = |, p = 2 in (b). In (b), we also show the first dipole profile 
fx(y) that, instead of ( 16. ip . satisfies the anti-symmetry conditions at the origin, 

(6.2) /( ) = f(0) = =^ = _/(„). 

Note an important feature of such compactly supported profiles that is seen in the fig- 
ures: by ( 12.71) . their mass must be zero. This necessary condition essentially "deforms" the 
VSS similarity profiles, so that it gets difficult to distinguish in Figure M their Sturmian- 
like properties on the numbers of dominant extrema and transversal zeros (if these apply 
at all). Note that the orthogonality property in (12. 7p is perfectly valid for the eigenfunc- 
tions (14. 3 p for n = (see (14.8 j) ) . which made it possible to develop the above branching 
theory. 
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6.2. p-bifurcation branches: numerics. Consider the semilinear case n — 0, which 
is known to be simpler, but correctly describes the expected general behaviour of the 
p-branches, at least, for sufficiently small n > (bearing in mind, that a continuous 
"homotopic" deformation as n — > is observed in a number of papers mentioned above). 
Thus, Figure [10] illustrates this subcritical case for even symmetry conditions (16. ip at 
the origin. The branches are seen to remain distinct, which contrasts markedly with the 
supercritical case [TO] [12] , where the branches are increasing with p, so that the P2 and 
higher branches "intersect" the vertical po branch, {p = 3} at points (profiles) / with the 
zero mass as in (12.71) . 

In Figure [TOT A) . we observe a strong, almost vertical, growth of these p-branches, which 
bifurcate, respectively, at 

(6.3) P3 = i + _^ = j j p 4 = i + _J_ = 2 f p 6 = i + _|_=£. 

This is not surprising, since the ODE (13. ip for iV = 1 and n = assumes, as p — > 1~ , 
balancing the terms 

(6.4) ... +(\f\P-if)»+ ... +^-L_ / = o =► f = Cf, where C7(p) ~ (p - 1)"^ 

(the scaled function f{y) is then supposed to be "almost" uniformly bounded, probably 
up to slower factors). Therefore, by (16. 4p . f(y) has a super-exponential growth as p — > 1 + . 
Since the bifurcation values in (16 .3p are already sufficiently close to 1 and the bifurcations 
are subcritical, these explain such a strong growth of all the p-bifurcation branches in 
Figure HDTA). 

In the nonlinear case n > 0, such an convincing justification of the general p-diagram 
is not available. Indeed, as we have shown in the previous section, in the simplest case 
/ = 0, i.e., p = po, the bifurcation of this vertical (in p) branch occurs from a nonlinear 
eigenf unction, which is the scaled source-type profile of the nonlinear thin film operator. 
Other nonlinear eigenf unctions of the thin film operator are still unknown possibly ex- 
cluding the second dipole-like eigenfunction. We expect that the discrete nature of the 
p-bifurcation branches discovered in [12] for n = remains valid for small n > 0, where 
the nonlinear branching points cannot be calculated explicitly as in (I2.10p and follows 
from a complicated nonlinear eigenvalue problem for the thin film operator. 

7. FBP: LOCAL BEHAVIOUR OF RADIAL SELF-SIMILAR PROFILES 

7.1. Interface conditions in the radial setting. As in [10], for the FBP, we need to 
look for profiles f{y), which vanish at finite y — yo > and describe asymptotics of the 
general solution satisfying the zero contact angle and zero- flux conditions: as y — > y$ , 

(7.1) /(y)->o, /'(v)-x), -^(-^(yN-ifyy + (\f\p-i f y ^o. 

The classification below also applies to the corresponding bundles occurring for the Cahn- 
Hilliard equation with n = 0, [121 § 2]. We assume that the conditions (I1.8P hold (as in 
the blow-up case [9]). 
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FIGURE 10. The bifurcation p-diagram and associated plots for the Cauchy 
case when N = l,n = 0. (A) shows the p-bifurcation branches emanating from 
the critical exponents p = pi = 1 + on the p-axis. The first four (even) 
branches I = 0, 2, 4, 6 are plotted. (B) illustrates the monotonicity of the mass 
of solutions in the critical case p = po = 3, whilst (C) shows selected profiles on 
the four branches in (A) that have ||/||oo = !■ 
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7.2. Two-dimensional asymptotic bundle of similarity profiles. The derivation 
of this bundle is standard and coincides with that in [TUT § 3.2]. Namely, as above, on 
integration, we obtain the ODE ( 13. 2D . Therefore, for n G [0, |), p > |, the two-parametric 
bundle of solutions is (cf. j5j [13] ) 

(7.2) /(y) = C (y - yf - A(y - y)*(l + o(l)), 

where yo > and Co > are arbitrary parameters and the correction term depends upon 
the value of iV and n, namely: 



a) for AT = 1, q = 5 - 2n, A 



= L 




q = 5 - 


9 = 


3. 


A = 


9 = 


3, 


A = 


9 = 


5 


- 2n, 



1-n 



(5-2ra)(4-2n)(3-2n) ' 



ifn<l, 

(b)foriV>2, < g = 3, A = - ^ if n = 1, 

, _ PvoC^- n) if n > 1 

^ ~~ (5-2n)(4-2n)(3-2n) 11 " ^ 

This expansion exists also for n = and has nothing to do with the CP exhibiting 
infinite propagation. It can be also used in the FBP posed for the Cahn-Hilliard equation 
with zero contact angle and zero-flux conditions. Proving such expansions demands a 
rather involved application of Banach's contraction principle; see also [5l [13]. 

8. FBP: SOURCE- TYPE SIMILARITY PATTERNS IN THE CRITICAL CASE p = p 

In the critical case p = p , the ODE (13. ip can be integrated once reducing the radial 
ODE to a third-order equation of the form 

(8.1) \f\ n {f" + /')' - (\f\ p °- l f)' ~ 0yf = 0, where = 

At the interface, we take conditions in ( 17. ip written as 
(8-2) /(y ) = /'(yo) = 0, 

and we complete the problem statement by taking the symmetry condition at the origin 
(13.111) . The ODE (18. ip itself then implies the second symmetry condition f"'(0) = 0. 

The mass M of / is a parameter, which is useful in distinguishing the solutions of (18.11) 
subject to (18.21) and (13. lip . Here, we set 

yo 

(8.3) M = fy N ^f(y)dy. 

o 

As another parameter, we can take /"(yo) or C = — |/"(yo), corresponding to the bundle 
( 17. 2p . There are also other choices of parameters, such as {/(0), /"(0)}. 

For the FBP, the local parameters {Co, yo} determine a two-parameter shooting problem 
in order to attain the single symmetry condition at the origin ( 13. lip . The global mass- 
parameter M is useful in classifying our solutions as bifurcations from critical mass values 
associated with non self-similar steady states. 

Thus the statement of the FBP comprises the ODE (18. ip with the symmetry condition 
(13. lip within the bundle (17. 2p with two parameters. Therefore, we expect a countable set 
of continuous families of solutions to exist, which can be parameterised relative to yo or 
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with respect to the mass M of the profiles. The asymptotic structure of the first (stable) 
branch of such similarity profiles is easier to detect. 

8.1. Continuous mass-branches of solutions of the FBP. For simplicity, we again 
consider the ID case. The global structure of similarity profiles does not essentially 
depend on N. In the critical case (jl.6p . i.e., p = n + 3 for N = 1, the ODE (18. ip takes 
the form (I3.10p . where at the interface point f(y) is assumed to belong to the bundle 
(17. 2p . Note that this equation is obviously non-variational. The existence of a continuous 
set of solutions with any sufficiently small mass is proved by a shooting argument exactly 
as in [12} § 5]. The only difference is that all the positive large solutions as y — > +oo are 
concentrated in a three-dimensional bundle around the profile 

f*{y) = y^V^ny), 

where the oscillatory component ip*(s) is a periodic function of changing sign of a certain 
autonomous ODE that is easy to derive. 

8.2. Numerical construction of the global similarity patterns. The boundary- 
value problem is now (18. ip written as 

(8.4) r + ^(r - \ /') -poiMf - (3yf\f\- n = o, 

together with (18 .2p . (I3.1ip . and (18. 3p . For fixed N and n, the numerical results presented 
below suggest that there is a countable number of solutions for a given y . We denote 
the profiles with positive mass as fk where the index k = 1,2, ... represents the number 
of sign changes (k — 1) of the profile over the interval [0, y ] (it also being related to the 
number of maxima and minima). The results for the base profiles f\(y) are presented in 
Figure [TT1 (similarity profiles) and Figure [121 (^/n-bifurcation diagram), whilst Figure [TBI 
gives the corresponding results for the next profiles f2(y) in the set. We remark that there 
are the corresponding reflected profiles —fk{y) which have negative masses. 

8.3. Asymptotic expansions. 

8.3.1. The small mass limit M — > 0. As in [TUJ § 4.2], we present an approach for the 
asymptotic expansion in the limit of small mass. This limit is easier for n = [T21 
Prop. 6.2], and, after rescaling, / converges to the fundamental similarity rescaled profile 
F defined in (14. ID as the first (with the eigenvalue Ao = 0) normalized eigenfunction of 
the linear non-self-adjoint operator B there. For n > 0, we identify the corresponding 
zero-mass limit as follows. Let us introduce the small parameter 

(8.5) e = M^ 
and perform the scaling 

2N nN nN 

(8.6) f (y) =£*+nN g(z), y = £2(4+niV) ^ y Q = £ 2(4+niV) Zq 
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FIGURE 11. Illustrative numerical solutions of the base / = f\(y) global sim- 
ilarity profiles in the cases N = 1,2. Parameter values of yo = 1 and 10 were 
used for the support length together with the selected values n = 0, 1, 1.5 for 
the index n. (A) and (B) show the one-dimensional case N = 1, (C) and (D) 
the two-dimensional case N = 2. In both dimensions, profiles when yo = 1 were 
not satisfactorily obtainable for smaller n values less than 1 (due to their values 
being below scheme tolerances). 



Then the problem (j8.ip ~ (j8.3p becomes 

i\9\ n (g" + ^9')'-^ = edging)' 
\ g (z )=g'(z ) = 0, g'(0) = 0, 1 = J c 



1.7) 



z o „N-1, 



- " g(z) dz, 



where ' now denotes 4-. We pose the regular expansions 



g = go -egi + 



ez + ... as e — > 0, 
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FIGURE 12. Parameter plots of the mass M = Mi of the base profiles solutions 
/ = fi(y) with yo for selected n values. (A) gives the one-dimensional case N = 1 
and (B) the two dimensional case N = 2. 
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FIGURE 13. Parameter plots of the second profile solutions / = f2(y) for 
selected n in one-dimension N = 1. (A) gives the profiles for yo = 10, whilst (B) 
gives the mass as the support yo varies. 



to obtain the leading-order problem 

to| n « + ^o)'-/^o = 0, z>0, 



(8.8) 



M4) = M) = 0, 9'o(0) = 0, 1 = z^goiz) dz. 
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The first-order problem, with z^ = — 9 }^ a l , is 



.9) 



g \ n {g" + ^g'l)' + (n - i)Pz 9l = -p g' \go\ po ~\ 

gi (z° ) = g[(z°) = 0, g[(0) = 0,0 = z N ^ gi {z) dz. 



These leading and first-order problems coincide with those in the unstable case de- 
scribed in [TUJ § 4.2.1]. As such the details will not be repeated, other than to correct 
typographical errors in (4.18) therein, which should read 

4i = - ^ im £ )iN+A)] * p»(&)> where = -iiW 12 , ft w = -u&r- 

We remark that the zero mass limit requires the support to vanish (as shown through the 
scaling (18.61) ). The corresponding limit of vanishing support but with the mass remaining 
non-zero does not possess non-trivial solutions in contrast to the unstable case. This is 
consistent with the physical interpretation of the terms on the RHS in (II. ip . since the 
effect of the second order term (the gravity) is no longer opposed by the fourth-order term 
(the surface tension) as in the unstable case (11. 2p . 

8.4. The large mass limit M — > oo. This limit occurs simultaneously with increase in 
support length y . We thus introduce the scalings 

7VH — 

(8.10) y = y z, f(y)=y^ 1 g(z), M = y ^m, 

so that (pHft - fjO]) becomes 

(8 n , , „ u 2(1+I ^ } (g" + ^g')' = Pzg\g\-* + Po \gfig', 




- g'{l) = 0, g'(0) = 0, m = z N - l g(z) dz, 

where ' again denotes 4-. This gives a singular perturbation problem in the limit yo — > oo, 
comprising an outer region < z < 1 together with an inner region near z = 1. 

(I) Outer problem. We pose the regular expansions 

g = g + o(l), m = m + o(l) as y -» oo, 
to obtain the leading-order outer problem 

I3zg \g \- n + Po\g \^g'o = 0, 

^o(l) = = 0, ^(0) = 0, m = z N - l g (z) dz. 
We thus obtain the explicit solution 

(8-13) go=[^{l-z-)]^, m Q = \{^]^ T 0^. 

V P q — 1 2 I 

This solution for g does not satisfy the condition g' (l) = and thus we require an inner 
region near z = 1. We note that this leading order outer solution is common to all profiles 
fk in this large mass limit. 
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(8.12) 



(I) Inner problem. Near z = 1 we introduce the scalings 



(8.14) z = l-5Z, g = 5^G, 



rajV+4 



where dominant balance in (18. lip gives the small parameter 5 = S(yo) = y nN+3 , and for 
Z = 0(1) we pose 

g = G + o(l) as y -)> oo, 
to obtain the leading-order inner problem 



(8.15) 



'G^ = -/3Go\G \- n +po\G \NG' , 
Go(0) = G'o(O) = 0, G ~[^-^]^ asZ^ 



where ' is now The last condition arises from the matching with the outer solution 
(18.131) . Consistent with the set fk, we anticipate a countable set of solutions denoted by 
Gqi-, k — 1, 2, ... , and distinguished by the number of sign changes (G fc having k — 1 sign 
changes). Numerical solutions for Go are shown in Figure [TH where the first two profiles 
are shown for the parameter value n = 1 in the N = 1 and N = 2 cases. 
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FIGURE 14. Numerical solution for the leading order inner problem in the 
large mass limit. Illustration of the first two profiles in a sequence of increasing 
changing sign profiles. The one-dimensional case N = l,n = 1 is shown in (A), 
whilst (B) shows the two-dimensional case iV = 2, re = 1. The curvature values 
at the origin may be used to distinguish the profiles and for comparison, we 
record the values G' \(0) = 0.2917, G" ' 2 (0) = -0.0960 for (A), whilst we have 
G' \(0) = 0.2130, G ' 2 (0) = -0.0667 in (B). 



2N 



9. FBP: ON COUNTABLE SETS OF p-BRANCHES OF SIMILARITY PATTERNS FOR p ^ p 

Similar to Sections H] and [5] for the Cauchy problem, we now intend to develop an 
analogous analytical approach to show existence of p-branches of similarity profiles in 
the FBP setting. This will demand rather unusual a special "spectral theory" for the 
corresponding linear problem for n — 0. 

9.1. The origin of countable p-branches. In view of essential differences and addi- 
tional difficulties arising for the FBP, we restrict ourselves to the simplest case N — 1, so 
that the FBP setting includes standard three free-boundary conditions 

(9-1) f(yo) = f\vo) = (\f\ n f"')(yo) =0 at an unknown boundary y = y > 0, 

accomplished with two symmetry (16.11) or anti-symmetry ones (16.21) at the origin. Overall, 
for the ID fourth-order ODE 

(9.2) A+(/) = -(\f\ n f"J + Pf'y + af + (l/r 1 /)" = 0, ye (0, y ), 

where a = 2p _(n+2) and (3 = 2 [2p-"n+2)] > these give five conditions plus an extra free 
parameter y (a nonlinear eigenvalue). This looks like a correctly posed problem, which, 
in the standard analytic setting, could not have more than a countable set of solutions, 
or a finite number of uniformly bounded ones. 

As usual, we next need to consider the corresponding ID pure TFE (I4.25p . for which, 
in the FBP setting there appears the following nonlinear eigenvalue problem (cf. (I4.26P ): 

rp ^ f-(l/l7'")' + 1 f i 2//' + a/ = on (y_, y+ ), 

1 } I / = / , = |/| n r = Oat y = y±. 

Note that, unlike the CP one (14. 26 p . the space of eigenvalues is three parametric, and, 
besides the usual nonlinear eigenvalue «6l includes two free boundary positions y±. In 
the two basic simpler cases with the symmetry (16.11) or antisymmetry (16. 2p . the eigenvalue 
space is 2D: 

(9.4) Eigenvalues: fi = (a, y ) T e M 2 , 
to which we concentrate upon in what follows. 

9.2. n = 0: first aspects of linear "Hermitian spectral theory". For n = 0, the 

nonlinear eigenvalue problem (19. 3 p becomes a linear one for the operator B in (14. ip : 

(9.5) = -V> (4) + | ip'y + \ if) = A^> (A = \ - a), ip = ip' = = at ?/ = t/ ± . 

Even in the simpler case (19 .4p . this is not a standard spectral problem, and, moreover, 
it is not clear whether this can be attributed to such classes. Let us comment that each 
"eigenfunction" ipk(y) is supposed to be defined on its own interval (yk~,yk+), but using 
eigenfunction subsets together with notions of completeness, closure, etc. may not make 
sense. 
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We now discuss some particular aspects of the problem ( 19. 5 p and restrict our analysis 
to the first eigenvalue Ao = and the corresponding even eigenfunctions. On integration, 
the ODE becomes of the third order, 

(9.6) -r' + |# = 0on (0,y„), V>(0) = 1, V'(0) = °5 V>(yo) = ^(Vo) = 0, 

where for convenience we fix the normalization ip(0) = 1 and take symmetry conditions 
at the origin (the subscript + then being dropped on the domain end point). We can 
prove the following: 

Proposition 9.1. The y^-eigenvalue problem {\9.6i) . corresponding to Ao = 0, admits a 
countable set of eigenfunctions {/0o > k > 1} defined on the intervals {(0, y$ )}, and, 

(9.7) as k -> oo, yj fc) ~ fct ->• +oo and 4> { Q k \y) ^ F(y), 

where F(y) = tp^iu) is the first eigenfunction o/ £/ie rescaled operator B m ( ^7 ) 
defined for the CP (i.e., in the whole R). 



Proof. A proper solvability of the problem ( 19 .6p follows from the known WKBJ-type 
asymptotics of solutions of this ODE for large y: 

(9.8) ^)~e^ 4/3 =► a 3 = 1(f) 3 . 

4 

This gives three roots: the real positive one ao = 3 • 4~3 and two complex: 

a±= (- |±i^)a . 

A full WKBJ expansion includes also a slow growing algebraic multiplying factor y^ 1 ^ 3 
(not important for the final estimates): as y — > oo, 

(9.9) if>(y) ~ y-ijde^ 473 + C 2 e"^ 4/3 cos (^yf + C 3 ) } + ... , 

where Ci i2i3 are real constants. Solving the problem at the free boundary y ^> 1, we see 
that acceptable roots are concentrated about the roots of the cos or sin functions, so that 
yf ] satisfy 

(9.10) cos(^(y«)i + C7s)«0 =► ^(^l^vrd + ^-Ca, 

whence the estimate in (19.71) . Since y^ — > oo, the convergence to the rescaled kernel F 
in the CP is obvious (see Figures below for a further justification). □ 

In Figure [T5| we show the first positive eigenfunction ^q 1 ^ of (19. 6p . Figure [TBI shows first 
four eigenfunctions of (19.61) . This illustrates the convergence to F in (19. 7p . where ip^\y) 

(5) 

is already very close to F, so that the next i/jq is difficult to detect numerically. 
For A 7^ 0, the corresponding statement to ( 19. 6 p is 



(9-H) 



y 

-ip"' + \y^ = \j if;(s) ds on (0, y+), 

#)) = i, ^(o) = o ; Tp{ y+ ) = Tp'{y+) = r , (y+) = o- 
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We distinguish the sets of even eigenfunctions {tp = ip m , X = X m ,y+ = Vm\i using the 
subscript m = 0,2,4,.... Within each set with m fixed, we have a countable set of 
eigenfunctions {ip m = ipm\y m = yln},k = 1,2,3,.... The first three eigenfunctions for 
the second m = 2 and the fourth m = 4 sets are shown in Figure [T71 



9.3. Nonlinear "Hermitian spectral theory". Consider now the nonlinear eigenvalue 
problem ( 19. 3 p for A = in the setting as in (19. 6p . 

(9.12) a = ^ -\f\ n f'" + ^yf = 0, f = /' = at the interfaces. 
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FBP,n=0: eigenfunctions\|/j, k) (y), k=1,2,3, A, =-0.354... 




12 14 



(a) four eigenfvmctions 
FBP, n=0: eigenfunctionsi|/^ k) (y), k=1,2,3, A. 4 =-0.617... 




(b) enlarged zeros 

FIGURE 17. First three profiles for the second and third even eigenfunctions 
m = 2,4 of f[97TT]t . 



Obviously, the first eigenfunction ■^q 1 '' is positive and is the classic one first obtained in 
[5] (see also [T3] for N > 1). However, we claim that, similarly to Proposition 19. 1[ the 
following holds: 
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FBP, n=1 : two eigenfunctions\|/jj , '(y), k=1 ,2, ^=0 






y<°°U.95.. 


Cauchy problem 









FIGURE 18. Two nonlinear eigenfunctions of <|9. 12 j> for n = 1. 



Conjecture [9J.2. Besides the positive solution [51 [13], in the oscillatory range n < 
given in {3.1), the problem <\9.12b admits a countable set of sign changing patterns f( k > = 
if>o (y), k = 1,2, ... , where each k-th one has precisely k zeros for y G (0, y^}- 

A rigorous proof is still incomplete since it requires detailed knowledge of oscillatory 
structures near interfaces such as ( 13. 5p . This should play a similar role to the linear 
expansion as in ( 19. 9ft . Such a deep understanding of the nonlinear expansion is still not 
achieved. 

In Figure [TBI we show the first two eigenfunctions of the problem ( I9.12p for n = 1. 
It is difficult to demonstrate more functions, since the next ones are already very close 
to the Cauchy profile denoted by ■ The oscillations of this CP-profile F(y) are 
presented in Figure [191 m between of the humps of which the interfaces of further nonlinear 
eigenfunction ij)^\y) are assumed to be situated. 

For A 7^ 0, the corresponding statement to ( 19 . 1 2[) is 



1+nA 



(9.13) 




yf = X I f( s ) ds ° n (o,y+), 

= 0; f(y + ) = f>(y + ) = \f\"f'"(y + )=0, 



with A 



l-q(n+4) 



. Again the sets of even eigenfunctions are denoted by {/ = ip m , A 



A m , y + = y m }, using the subscript m = 0, 2, 4, Within each set (i.e. m fixed), we have 

a countable set of eigenfunctions |-0 m = ipm\ y m = yffl^ ,k = 1, 2, 3, ... . The first three 

eigenfunctions for the second m = 2 and the fourth m = 4 even sets are shown in Figure 
[201 for the case n = 1. In the nonlinear case n > 0, unlike the linear one n = above, 
convergence of numerical methods become much more slower. 
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9.4. Comments on n-branching and p-bifurcations. These properties are assumed 
to be similar to those developed earlier for the CP-setting. However, there are essential 
difficulties even in doing some formal computations. 

Overall, we expect that the nonlinear eigenvalue problem f lQ. 1 2[) (see Conjecture 02) 
admits n-branching as n — > + from linear eigenf unctions from Proposition 19.11 As usual 
"variations" of the free boundary then must be taken into accound, which does not lead 
to any extra difficulty. 

Furthermore, the general VSS problem (19. 2 p possesses p-bifurcation branches obtained 
via nonlinear bifurcations, whose theory is developed along the same lines used in Section 
15.21 In particular, as in (15. 9p . the first bifurcation exponent will be 

Similar to the analysis in Section I5.2[ further study is necessary to check whether a 
"nonlinear bifurcation" occurs at p = po (possibly not as for n = 0), so further critical 
exponents a = an(n) should be revealed. 

We also expect that the above set of nonlinear eigenfunction-eigenvalue {fi(y), an(n)} 
pairs for the CP for the TFE (ll.9p also play a role for the FBP. More precisely, it is 
expected that, for any given CP profile fi(y), with a = ai(n), there exists a countable 
set of FBP profiles {fi(y), ai(n)}, which are defined on the expanding intervals [0,yj), 
where y J — > y^ (the interface of the CP-profile) as j — > oo. Eventually, there holds: 

(9.14) fi(v)->fi(v) as j^oo 

uniformly on compact subsets. In other words, each CP profile fi(y) can be arbitrarily 
closely approximated by FBP ones with a finite number of sign changes. The convergence 
(I9.14p is a difficult open problem, both analytically and numerically even for the first 
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FBP, n=1: eigenfunctions\|/^ ',k=1,2,3, A, =-0.19... 
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(a) four eigenfvmctions 
FBP, n=1: eigenfunctionsi(/^ ) ,k=1 ,2,3, A, =-0.31... 




-1 ' 1 1 1 1 1 1 1 1 1 1 

0123456789 10 

y 

(b) enlarged zeros 

FIGURE 20. The first three profiles for the second and third even eigenfunctions 
m = 2, 4 of (|9.13p for n = 1. These numerical solutions of (|9.13p used the 
regularisation (13. 12ft with 5 = 0.1. 



values of I. Of course, (I9.14p is associated with the discovered oscillatory properties of all 
the VSS CP profiles. 
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Again, we mention that the FBP problems have a different nature and are more difficult 
than the CP, since a new free parameter the interface location y Q > (an extra "nonlinear 
eigenvalue") appears in the mathematical setting. 

10. Discussion 

We began in Section [2] with the similarity analysis of global, source-type (for p = po) 
and very singular (for p ^ po) solutions of the TFE with a stable parabolic term (11. ip . 
Treating the Cauchy problem (respectively, the FBP), we first described in Section [3] for 
the CP (respectively, [7] for the FBP) the local asymptotic properties of solutions near 
interfaces. While the FBP asymptotics turned out to be standard and reasonably well 
known since the 1990s, it is important that, for the CP, we detected the nonlinearity range 
n G (0, rth) (rih is the point of a heteroclinic bifurcation for a nonlinear related ODE), in 
which the rescaled ODE exhibits a unique stable periodic motion describing, as expected, 
generic changing sign properties of more general solutions. 

In Section [3j we studied the CP in the critical case 

p = p = n + l + ^, 

where we detected continuous branches of similarity profiles. In Section HI we developed 
n-branching theory of similarity solutions of the ID pure TFE 

u t = -(\u\ n u xxx ) X} 

where we described branching of nonlinear similarity profiles from eigenfunctions of a 
linear rescaled operator at n = 0. This allowed us in Section [5] to reveal a countable se- 
quence of critical exponents {pi} of the original stable TFE (II. ip and to describe similarity 
solutions for p ^ p . 

After a detailed study of the CP, we returned to the FBP setting. We studied in 
Section [S] various branches of similarity patterns for the FBP in the critical case p = p 
and extended some of the results to p ^ p in Section O 

By comparing the similarity patterns of the CP and the FBP, a striking "limit" property 
emerges: the infinitely oscillatory patterns of the CP are the limits of FBP-patterns with 
a finite number of sign changes. Naturally, this is required to take into account sign 
changing patterns of the FBP, which have not been previously studied in any detail. In 
a sense, the above limit property can be considered as a certain definition of solutions of 
the CP (besides the already existing ones via maximal regularity at the interfaces or via a 
smooth analytic "homotopy" deformation to the bi-harmonic equation u t = —A 2 u, which 
turned out to be a good approximation of the TFE for small n > 0; [10J). 

Thus, the goal of the paper was to describe some leading key ideas concerning (i) 
formation of similarity solutions for the stable TFEs and (ii) extra relations between 
the CP and the FBP. Some of the most difficult conclusions remain formal, the related 
mathematics turns out to be very difficult, and we have posed several open problems for 
future research. 
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